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ABSTRACT 

Based upon a novel paradigm of the galaxy formation under UV background, the evo- 
lutionary bifurcation of pregalactic clouds is confronted with observations on elliptical 
and spiral galaxies. The theory predicts that the dichotomy between the dissipational 
and dissipationlcss galaxy formation stems from the degree of self-shielding from the 
UV background radiation. This is demonstrated on a bifurcation diagram of collapse 
epochs versus masses of pregalactic clouds. Using the observed properties, the collapse 
epochs are assessed for each type of galaxies with attentive mass estimation. By the 
direct comparison of the theory with the observations, it turns out that the theoreti- 
cal bifurcation branch successfully discriminates between elliptical and spiral galaxies. 
This suggests that the UV background radiation could play a profound role for the 
differentiation of the galactic morphology into the Hubble sequence. 



Key words: galaxies: formation — radiative transfer 



molecular processes 



1 INTRODUCTION 
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' ' I ( |Silk 1977D ). The theory predicts that the galac- 
tic scales are basically determined by atomic cooling of hy- 
drogen and partially helium. The pregalactic evolution is 
elegantly summarized on the cooling diagram for virialized 
objects. Moreover, the origin of the Hubble type has been at- 
tributed to the dissipativeness of the collapse, which is regu- 
lated by the efficiency of star formation. If the star formation 
proceeds after most of the gravitational energy is dissipated, 
the pregalactic clouds will evolve into spiral galaxies. This 
is a p aradigm of th e so-called dissip ati onal galaxy formati on 
(e.g. flLarson 1976] ); ( |Carlberg 198^ ); ( |Katz fc Gunn 199l| )). 
On the other hand, an early star formati on episode leads to 
the d iss ipationless galaxy format ion (e.g. ( Aarseth & Binney 
197§); flAguilar fc Merritt 199G-)), ending up with the for- 



berg 1981|); (palla, Salpeter fc Stahler 1983|); flSusa, Uehara 



mation of elliptical galaxies. However, the key mechanism 
which physically controls the star formation efficiency has 
been hitherto unsolved. 

Very recently, Susa & Umemura (2000) (hereafter SU) 
propose that the self- shielding against UV background ra- 
diation regulates the star formation in pregalactic clouds. 
The star formation processes in primordial gas have been 
explored by many authors with ab initio calculations (e.g. 



* e-mail:susa@rccp. tsukuba. ac.jp 

t e- mail : umemuraOrccp . tsukuba. ac.jp 



flNishi et al. 1998Q ; ( |Nakamura fc Umemura 1999[ )). The 
key physics is the radiative cooling by H2 line emission, be- 
cause H2 is the only coolant for primordial gas at T 10 4 
K. SU have studied the efficiency of H2 cooling in collapsing 
clouds exposed to UV background radiation, because it is 
significant after the reionization of the univers e, probably 
at 2: ^ 10 (( Nakamoto, Umemura, fc Susa 1999] ), and refer- 
ences therein). SU have found that if a cloud undergoes the 
first sheet collapse at higher redshifts (z <; 4), then the cloud 
is quickly shielded against the UV background and conse- 
quently cools down due to the efficient formation of H2 . Re- 
sultantly, it leads to an early burst of star formation. On the 
other hand, the shielding is retarded for a later collapsing 
cloud at z S5> 4, resulting in the dissipational galaxy forma- 
tion. As a result, the bifurcation branch of the self-shielding 
is corresponding to the boundary between the dissipationless 
and dissipational galaxy formation. In this Letter, this novel 
bifurcation theory is confronted with the observations of el- 
liptical and spiral galaxies to elucidate whether the theory 
is practically successful or not. 



2 BIFURCATION THEORY 

In Figure 1, we show the bifurcation theory. Here, the cos- 
mological parameters are assumed to be = 0.3, h — 0.7, 
and Q\,h 2 = 0.02 with usual meanings. Originally, the sheet 



© 2000 RAS 



2 Susa & Umemura 



collapse was pursued with two initial parameters, i.e., the 
mean density (nini) and the thickness (A). Here, the ini- 
tial parameter space is translated into the baryonic mass 
[Mb = (47r/3)ni n i(A/2) 3 ] and the collapse epoch (z c ) by as- 
suming the initial stage is close to the maximum expansion 
of a de nsity fluctuation. In the region (a) in Figure 1, a pre- 



tion of stellar component in the total baryonic mass. The 
mass-to- light ratios M,/L are obtained theoretically as well 
a s observationally. Based on the chemical evolution theory 
(( |Kodama fc Arimoto 1997| )), the M,/L at B-band for el- 
lipticals is 4 — 9. This valu e is consistent with t he observa- 
tional estimates of 4 — 7 ((Bertola et al. 1993] ); (Pizzella et 



galactic cloud is self-shielded against the external UV in the al. 1997)). In this paper, we adopt the luminosity- dependent 



course of the sheet collapse, so that the cloud cools down 
below 10 3 K and undergoes efficient star formation. Hence, 
it is expected to evolve into an early type galaxy with a 
large bulge-to-disk ratio (B/D) due to the dissipationless 
virialization. In the region (b) , the cloud is not self-shielded 
during the sheet collapse, but will be self-shielded through 
the shrink to the rotation barrier. This leads to the retarded 
star formation, and thus the virialization would proceed in 
a fairly dissipative fashion. As a result, a late type (small 
B/D) galaxy would be preferentially born. The region (c) 
represents the forbidden region of the collapse due to the 
Jeans stability. 

SU have considered only baryonic component, because, 
in the later collapsing phase, the cooling sheet is dominated 
by baryons, not by the diffuse dark matter. Then, if the UV 
background is constant, the bifurcation mass scale is given 



Msg* = 2.2 x IO^Mq 



(1) 



where I21 is the UV background intensity in units of 
10 _21 erg s _1 cm _2 str _1 Hz _1 . However, the dark matter po- 
tential may affect the evolution in the early phase of col- 
lapse. If we take this effect into accoun t max imally, the 
Jeans length is reduced by a factor of yQ^jQ.. With this 
Jeans scale, the bifurcation mass is changed to be M§^ n = 
(fib/fi) 1 ' 2 Mg^ x . Thus, the practical bifurcation mass would 
be between M^ ax and Mgg n . 

Also, SU have assumed the UV intensity to be indepen- 
dent of time. Practically, the intensity seems to evolve. Here 
we include the effect of the evolution of UV background ra- 
diation. We assume J21 = 0.5 [(1 + z) /3] 3 for z < 2 and 
I21 = 0.5 for 2 < z < 4. This dependen ce is consisten t 
with the UV intensit y in the present epoch ((Maloney 1993); 
(Dove & Shull 1994)), and the valu e inferred from the QSO 



proxim ity e ffects at high redshifts ((Bajtlik, Duncan, & Os- 
triker l]98§); Giallongo et al. 1996). As for z > 4, two ex- 



treme models are employed. The first one is (A) the exponen- 
tially damping model, I21 = 0.5 exp [3 (4 — z)] ((Umemura et 
al. 2000J)), and the second one is (B) the constant extrapola- 
tion model, /2i = 0.5. In Figure [D, these two extreme models 
are shown. The difference emerges especially at z c <; 7. 



ESTIMATION FOR GALACTIC MASSES 
AND COLLAPSE EPOCHS 



Here, b hacd upon the obacrvational data, we attempt to aa 

sess the total baryonic masses of ellipticals as well as spirals 
and their collapse epochs. Then, the observed galaxies are 
compared with the theory in the bifurcation diagram. 

To begin with, we evaluate the baryonic masses from B- 
band luminosities for ellipticals and from I-band luminosities 
for spirals with the mass-to-light ratio; Mb = [M*/L] Lj /», 
where M* is the total stellar mass and /» is the mass frac- 



values given by Kodama & Arimoto (1997). The M»/L at I- 
band for spirals is 3.37(Sa), 2.91(Sb), 1.79(Sc), and 1.33(Sd) 
based on the code developed by Kodama & Arimoto (1997) 
with the SI model in Arimoto, Yoshii, & Takahara (1992). 
These are also consistent with the observed data in Rubin 
et al.(1985). /, is theoretically calculated by the population 
synthesis model to be 0.963(Sa), 0.908(Sb), 0.462(Sc), and 
0.125(Sd), while basically /» = 1 for ellipticals. With the 
estimated baryonic masses, we can assess the total masses 
as M to t = M b (ft/fib). 

Next, we estimate the collapse epochs with the help of 
virial theorem. For the purpose, we use the observed 1-D 
velocity dispersion. We assume that the system is spher- 
ical and the dark halo has isothermal distributions after 
the virialization. We suppose a density perturbation as 
8%{r) — Sig(r), where r is the comoving radial coordinate 
and g(r) is a function which satisfies the normalization as 
jjs g{r)r 2 dr = 1, with 7? being the comoving radius of 
the perturbed region. 5i is the spatially averaged 5i(r) in the 
volume r < R. Summing up the initial kinetic energy and 
gravitational energy, we have the initial total energy of the 
perturbed region as 



5R 



{l + z 



)&(l 



<t> )F, 



(2) 



where z; is the initial redshift, <j> is the contribution from 
the peculiar velocities, and M to t is the total mass enclosed 
within the comoving radius 7?. The factor F is defined as 



/ -^5 / r 4 g{r),ir. 



(3) 



where g(r) = ^ r' 2 g(r')dr' '. In the linear regime, g{r) oc 
r V(n+3) around a density peak (( [Hoffman fc Shaham 198^ )), 
where n denotes the index of the CDM power spectrum. 
Then equation ([|) is readily integrated to give F = 5/(2— n). 
We assume n = — 1 throughout this paper, because it is the 
case for galactic scales in CDM cosmologies. The total en- 
ergy of an observed galaxy is derived in terms of the virial 
theorem; E h B = — |M to tcr 2 D , where am is the 1-D inter- 
nal velocity dispersion of the galaxy. Assuming the energy 
conservation, i.e., 7?i nt = 7? b s , we have 



(1 + ^)^(1 + 20/3) 
R 



5er 2 D _ F -i 



2GM U 



(4) 



On the other hand, the time when the outer b oundary 
of a p erturbation collapses is analytically obtained ((Peeblei 
1993; )) as 



t = 



R 



(GM tot ) 1/2 \2(l + z t )5 l (l + 2<t>/3) 



3/2 



(5) 



Finally, using equations (U|) and (H), we have the collapse 
epoch as t = irGMtotF 1 /5 ' cr 1D . By equating this time 
with the Hubble time tn(z c ), t can be translated into the 
collapse redshift, z c . Thus, it turns out that the velocity 
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dispersion is a key quantity to determine z c . We use the ob- 
served stellar velocity dispersion as <tid for elliptical galax- 
ies. Also, for some luminous ellipticals, am is est imated by 
the X-ray temperature Tx as axo = \J kTx / p,m v . The esti- 
mation by X-ray gives typically twice the optical estimation, 
being probably the maximal assessment of a 2 D . For spiral 
galaxies, we interpret the asymptotic rotational velocity into 
the velocity dispersion by the relation crxa = v TO t / v2. 

The basic data on the luminosities and the velocity dis- 
persions are taken from the table of Faber et al. (1989) for 
332 ellipticals. The galaxies are selected from the original 
sample on the condition that the surface brightness, angular 
size, and velocity dispersion are measured. The X-ray data 
for 12 ellipticals are adopted from the table in Matsumoto et 
al. (1997). Further, for 468 spiral galaxies, we combine the 
Fband luminosities from Mathewson et al. (1996) and the 
asymptotic rotation velocities from Persic & Salucci (1995). 
We pick up the galaxies whose asymptotic velocities [V aa as 
in Persic & Salucci (1995)] are given. Thus the number of 
galaxies is smaller than the original one. 



4 THEORY VERSUS OBSERVATIONS 

In Figure ^| the bifurcation theory is confronted with obser- 
vations. In this figure, cosmological parameters are assumed 
as Q = 0.3, Ov = 0, h = 0.7, and Q. h h 2 = 0.02. All of these 
values are plausible in the light of the recent observations, 
although the value of A parameter is still controversial. The 
small filled triangles and open circles represent respectively 
the observed elliptical and spiral galaxies. The open stars are 
the X-ray luminous elliptical galaxies. We find that the two 
types o f galaxies are successfully divided bv the theoretical 



bifurca tion mass scale. This implies that the self-shielding 
against UV background radiation is practically related to 
the origin of the galactic morphology. The difference be- 
tween model A and B is not so significant, although model 
A predicts a bit more ellipticals than observed. Also, it is 
worth noting that few spiral galaxies reside above the bi- 
furcation mass scale, whereas there are a noticeable number 
of elliptical galaxies which seem to have formed at lower 
redshift epochs than predicted by the present theory. Since 
a condition for constant velocity dispersion gives a relation 
M oc (1 + z)~ 3 ^ 2 in this diagram, most of discrepant low- 
redshift ellipticals turn out to be relatively luminous galaxies 
with higher stellar velocity dispersions. Intriguingly, Gonza- 
lez (1993) has reported, using strengths of Hp and [MgFe] in- 
dices, the evidence for intermediate-age populations in ellip- 
tical galaxies. In addition, it has been argued that the galaxy 
merger with burst-like sta r formation can lead to the forma- 
tion of elliptical galaxie s ((Larson fc Tinsley 1976); B arnes & 
Hernquist 1991, 1996; flKauffmann fc Char lot "l998[ ) ) . Also, 
Shier & Fischer (1998) suggest, by studying stellar kinemat- 
ics of starbursting infrared-luminous galaxies with obvious 
morphological signatures of merger, that they can be pro- 
genitors of elliptical galaxies with high stellar velocity dis- 
persions. Thus, the discrepant low-redshift ellipticals might 
not be of primordial collapse, but could be a category of 
merger remnants. 

In Figure |^, la, 2a, and 3<r CDM density perturba - 
t ions normalized by COBE data ((Hu & Sugiyama 1996); 
(Bunn fc White 1997)) are also plotted. We find that ellip- 



tical galaxies form typically from ~ 4<r perturbations and 
spirals do from ~ 2a peaks. If we change the cosmology, 
the results alter to some degree, because the bifurcation 
mass, the CDM spectra and the data points of elliptical, 
and spiral galaxies are relatively changed. The main change 
of Figure comes from changing Q. The collapse epoch z c 
for the observational data is dependent upon Q as roughly 
(1 + z c ) oc Q-°' 74 for fi A = universe. This dependence 
is mainly due to the evaluation of dark mass in a galaxy. 
For instance, if we employ £1=1 (although it is unlikely), 
1 + z c becomes almost 2.5 times smaller than the estimates 
in Figure || In this case, the theory and the observational 
data seem perceptibly discrepant with each other. On the 
other hand, Figure |^ is insensitive to the cosmological con- 
stant parameter £Ia- In fact, even if we employ £1 = 0.3 
and £Ia = 0.7, the relative position between the bifurcation 
mass and observations remains almost unchanged, although 
the correspondence between CDM fluctuations and observa- 
tions shifts from ~ 4a to ~ 3a for ellipticals. In this case, the 
theory is in a good agreement with the observations. Fur- 
ther details of the dependences on the cosmological model 
will be discussed in a forthcoming paper. 

The present results are also intriguing from a view point 
of the statistics of galaxies. Bardeen et al. (1986) have shown 
that higher a peaks reside preferentially in denser regions 
rather than in low-dense regions. As a result, the so-called 
density morphology relation can be explained as a natu- 
ral consequence of the bifurcation theory. Furthermore, it 
is known that the specific angular momenta J/M of spi- 
rals are systematically greater by a facto r of <; 3 t han J/M 
of ellipticals in the same mass range ((Fall 1983)). Based 



upon the tidal origin of the angular mome ntum, we expect 
J/M oc (1 + z c )~ 1/5 V _1 , where v = 5/a (( [Heavens fc Pea- 
cock 198£ )). If ellipticals form from ~ 4a and spirals do from 

3(J/M) cll 

iptical - 



2a, then we anticipate ( J/M) sp - 1T 
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1 10 1 10 

l+z 1+z = 

(a) (b) 

Figure 1. The bifurcation diagram. The prediction of galactic evolution is illustrated on the plane of the baryonic mass (M^,) versus 
collapse redshift (z c ). Fig.[l]a corresponds to the UV evolution of model A (see text), and Figjljs for model B (see also text). In both 
of the figures, the half-tone shaded area (a) denotes the region in which pregalactic clouds evolve into early type galaxies. The shaded 
area (b) leads to the formation of late-type galaxies. The unshaded area (c) represents the region in which the density perturbations 
cannot grow into bound objects, due to the thermal pressure of ionized gas. Here we employ as the bifurcation mass, although 

the effect of dark matter may reduce the mass by factor (f2 b /f2) 12 at most. M shic i d denotes the mass above which the cloud is initially 
self-shielded. 
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Figure 2. Comparison of observed galaxies with the bifurcation theory. Small open circles represent spiral galaxies in Persic & Salucci 
(1995). Small solid triangles denote the E and E-SO galaxies in Faber ct al. (1989). Open stars are elliptical galaxies from X-ray 
observations in Matsumoto et al. (1997). The thick solid and dashed lines represent the bifurcation mass M^ x , respectively for model 
A and B, where both are identical at z c < 4. Three thin solid lines marked as la, 2a, and 3a are the prediction in the CDM cosmology, 
where a is the variance of CDM perturbations. The downarrow in the upper-left corner of the panel shows the maximal shift of the 
bifurcation mass from to M^ n when including dark matter. 
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